Vector interaction strength in Polyakov-Nambu— Jona-Lasinio models from 

hadron-quark phase diagrams 
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We estimate the vector interaction strength of the Polyakov-Nambu-Jona-Lasinio (PNJL) 
parametrizations, assuming that its transition curves should be as close as possible of the re- 
cently studied RMF-PNJL hadron-quark phase diagrams. Such diagrams are obtained matching 
relativistic mean-field hadronic models, and the PNJL quark ones. By using this method we 
found for the magnitude of the vector interaction, often treated as a free parameter, a range 
of 7.66 GeV" 2 < Gv < 16.13 GeV" 2 , or equivalently, 1.52 < Gv/G s < 3.2, with G s being the 
scalar coupling constant of the model. These values are compatible but restricts the range of 
4 GeV~ 2 < Gv i; 19 GeV -2 , recently obtained from lattice QCD data through a different mean- 
field model approach. 

PACS numbers: 12.38.Mh,25.75.Nq 

The hadron-quark phase transition is still a challenging 
task for both theoretical and experimental fields. From 
the theoretical point of view, the strongly interacting 
matter is treated by QCD. However, in the regime of 
low energies, QCD is nonperturbative and still difficult 
to solve for intermediate temperatures and chemical po- 
tentials, although lattice methods have been faced a huge 
progress in the last years [l[ . For such regime, it is useful 
to use effective models in the description of the quark 
matter that share the same features of QCD. For exam- 
ple, the MIT bag model [3, the NJL one 0, H and its 
version coupled with the Polyakov loop, named as PNJL 
model [f|. 

In this context such models are used to construct 
the QCD phase diagram |(|, where the different regions 
are identified as those in which the chiral symmetry is 
broken or restored. Studies in this direction were per- 
formed, for instance, for the linear a model 0, and the 
NJL one Q. The additional information about the con- 
fined/deconfincd phases can also be taken into account 
when the PNJL model is used to construct the quark 
phase transitions 0-HH , or even when the Polyakov loop 
is linked with the linear a model [ijj . 

In a very recent study [i~3| , it was performed a compar- 
ison among the hadron-quark phase diagrams generated 
by PNJL models, and those constructed by matching 
a large class of relativistic mean-field (RMF) hadronic 
models with four different parametrizations of the PNJL 
quark model (other studies based on this treatment can 
be found in Refs. fl4l - [il| ). The results shown pointed 
out to a difference between the phase transition curves 
due to the repulsive interaction of the RMF models. 

Based on these results, we propose the construction 
of a PNJL model that minimizes this difference. This 
will be done by including a vector interaction in the orig- 
inal PNJL structure. We name hereafter the resulting 
model as PNJLv model, and the strength of its vector 
interaction will then estimated in order to approximate 



the PNJLv transition curves to the RMF-PNJL ones as 
much as possible. 

Actually, the inclusion of vector interactions in effec- 
tive quark models was already discussed in the literature, 
see Ref. [|[ for a study in the context of the NJL model. 
It is known that the effect of the increase of the repul- 
sive interaction strength in the quark matter phase dia- 
gram is to shrink the first-order transition region as the 
chemical potential [i increases. Also the critical end-point 
of the transition moves to larger [i and lower tempera- 
ture T. Such effects, and further studies are reported for 
both NJL model [11 HI and the PNJLv one [lo|, Q1 Sl|- 
[23j | . where the vector interaction strength is often used 
as a free parameter. Therefore, taking into account the 
aforementioned effects, we furnish here a method to esti- 
mate this interaction strength, based on the RMF-PNJL 
hadron-quark phase diagrams. We also compare our re- 
sults with the recent ones, based on a mean-field calcula- 
tion related to QCD lattice data, proposed in Ref. [24| . 

Our starting point is the Lagrangian density of the 
two-flavor PNJLv model, that reads 



(1) 



with _D M = — iA^ being the covariant derivative, 
the gluon field, and m the current quark mass. The 
vector interaction strength is regulated by the param- 
eter Gy, and the Polyakov loop potential is given by 
U(T, fj,, $, $*). As in Ref. [l]|, we use four different ver- 
sions of this potential, namely_RRW06 HSU, RTW05 
[H, FUKU08 and DS10 g7|. 



From this Lagrangian density one obtains, following 
the procedure used in Ref. [26j | and taking into account 
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the new vector term, the grand thermodynamic potential, 
Q = U(T,/j,, $,$*) + G s p 2 s - ^ / E(k,M)k 2 dk 
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in the isospin symmetric system, in which E(k,M) = 
(k 2 + M 2 ) 1 / 2 , and with the degeneracy factor given 
by 7 = 12. The constituent quark mass is M = 
m — 2G s p s , and the quark density is obtained from 
p = —dil/dp, where p is related with Gy and p through 
ft = p — 2Gyp. Note that, the quark density can be 
also found by requiring that dVt/dp = 0. The func- 
tions T± = T±{k,T, ft, <&*) are defined by T± = 
F(k,T,fL,$,$*)±F(k,T,il,$,$*). 

The strength of the scalar interaction (G s ), the vacuum 
integral cutoff (A), and m are also parameters of the 
model. Here they are given by G s = 5.04 GeV~ 2 , A = 
651 MeV, and m = 5.5 MeV 

The inclusion of the confinement information, via the 
Polyakov loop, affects the statistical distributions of the 
PNJLv model in such way that the new Fermi-Dirac func- 
tions are now given by 



F(k,T,ft, *,$*) = 
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and F(k,T,fi,<$>,$*) = F(k, T, —ft, $*, $), with 
x = (E- ft)/T. 

The quark condensate ('0V') = Ps, and the Polyakov 
loop <!> are found by requiring that dCl/dp„ = 90/9$ = 
0, in the lowest order approximation [9j, |28|, that leads 
to <E> = <J>* . The explicit form of the equations of motion 
(EOM) of the PNJLv model, found by minimizing fi in 
respect to p, p s and <£, are: 



and 

where 
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and g2 = gi(k,T, —ft, $). Note that these EOM are the 
same PNJL ones, with the chemical potential p shifted 
by the vector interaction, as given by p. In the RMF 
case, such a shift also occurs in the value of the baryonic 
chemical potential due to the contribution of the mean- 
field value of the isoscalar vector meson field uj. 



It is important to remark that we are dealing with the 
simplest PNJLv version, in which color quark conden- 
sates are not being taken into account. Here, we are not 
considering any possibility of emergence of two-flavor su- 
perconducting color (2SC), or color flavor locked (CFL) 
phases [2i|[3c|. 

To construct the PNJLv transition curves varying 
the Gy parameter, we follow the procedure adopted by 
Fukushima [Io[ that uses the magnitude of the order pa- 
rameters p s and <I>, to define the transition temperature 
for each fixed p. In that work, the author uses the con- 
dition of p s /pT c = 1/2 to construct the phase diagrams. 
The transition temperature at p = found in that case 
is around T c (p = 0) = 200 MeV. Noteworthy to observe 
that in the construction of the quark phase diagrams, the 
approach using fixed values for p s is qualitatively equiv- 
alent, at least at moderated p values and for Gy ^ 0, 
to those based on the local maximum of dp s /dT and 
9$/9T, used in Ref. (2(|. In the Gy = case, the agree- 
ment between the crossover transition curves obtained 
with the two criteria is still better. 

Here, we use different values of p s / p™ c for each 
Polyakov potential in order to obtain a better agree- 
ment of T c (p = 0) with the lattice QCD results, 
T c (p = 0) = 173 ± 8 MeV [p. The adopted values 
arc Ps/pT" = 0.73,0.70,0.72,0.71, respectively for the 
RRW06, RTW05, FUKU08, and DS10 parametrizations. 
Furthermore, we still maintain the rescaling of the origi- 
nal parameters T and b of U{T, to T = 190 
MeV (RRW06, RTW05 and DS10), and b = 0.007A 3 
(FUKU08) also used in Ref. [lj]. 

The PNJLv diagrams constructed from the adopted 
method are displayed in Figs. H^JJJi. 

From Fig. [TJ one can see that there is a range of val- 
ues for Gy, at least in a certain temperature region, that 
makes the transitions constructed from the PNJLv model 
very close to those obtained via the RMF-PNJL match- 
ings, represented by the gray bands. For these match- 
ings, we have used a large class of RMF hadronic models 
coupled to the PNJL ones in which Gy = [l3j . 

The overlap between the hadron-quark phase tran- 
sitions provided by the PNJLv and RMF-PNJL 
models, is found for G V IN < Gy 
with G v m /G s = 1.60,1.52,1.60 
G™ AX /G S = 3.20,2.74,2.86,2.98, 
for the RRW06, RTW05, FUKU08, 
Polyakov potentials. This give us a total range of 
7.66 GcV" 2 < Gy < 16.13 GcV~ 2 . We can also define a 
temperature region of better overlap between the PNJLv 
results, and the RMF-PNJL bands, given by T < T MAX . 
The maximum temperatures are T MAX « 90, 80, 50, 80 
MeV, for the same aforementioned PNJLv models. 
Therefore, our findings can be useful for instance, in the 
study of proto-neutron stars (PNS) that are described 
at T < 50 MeV. Applications of the RMF-PNJL models 
to compact stars, have been done recently for the PNS 
evolution [321 . for quark [33[, and hybrid stars [lij]. 
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FIG. 1: (Color online) PNJLv phase diagrams results for different Gv values compared to RMF-PNJL models. The crosses 
indicate the critical end-points. 



In terms of the ratio Gy/G s , the total range ob- 
tained here is 1.52 < Gv/G s < 3.2. This result sug- 
gests that the Gy values calculated from hadron-quark 
phase transition curves, are in fact greater than some 
used in the literature, namely, Gy/G s = 0.25 J3f| and 
Gv/G s = 0.5 [36j |. In Rcf. I22j it was argued that the 
range of acceptable values in the non-local PNJLv model 
are 0.25 < Gy/G s < 0.5. This range is substantiated 
by a Fierz transformation of an effective one-gluon ex- 
change interaction, with Gy depending on the strength 
of the {/a(1) anomaly in the two-flavor model. Our val- 
ues are considerable above 0.5, which should not be un- 
expected, as the vector term in the PNJLv model, in 
our fitting procedure, mimics the repulsive short-range 
part of the nucleon-nucleon interaction present in the 
RMF model. The nature of the repulsive nuclear force 
is beyond one gluon exchange and acts between color- 
less hadron degrees of freedom, that essentially are of 
non-perturbativc origin, explaining the difference in the 
values of Gy. Consistently, our values are larger than, 
for instance, Gy/G s = 1 used for the PNJLv model of 
Ref. 0. 

It is important to stress that our estimated 
range for the vector coupling strength, is compati- 
ble with those found recently in Ref. |24j |. given by 
4 GeV~ 2 < Gy < 19 GcV~ 2 . Actually that range en- 
compasses ours. In that work, the authors estimated 
such range through a completely different way. They 
used the lattice QCD data of the diagonal and off- 
diagonal quark susceptibilities, as input in an effective 
QCD mean-field model approach, described by a La- 



grangian density of a two-flavor quark system interact- 
ing only via massive vector fields. By doing so, they 
assume that the vector part of the QCD interaction can 
be isolated, and treated in that approximation. Actu- 
ally, the authors considered a temperature dependent 
vector coupling, Gy = Gy(T). They obtained the range 
4 GeV" 2 < G V {T C ) < 19 GeV" 2 by assuming a transi- 
tion temperature of T c (p, = 0) = 170 MeV. In our calcu- 
lations we consider a fixed vector interaction strength. 

For the sake of completeness we also present in Fig. [2] 
the order parameters p s and and the quark density, 
obtained sclf-consistently from the PNJLv EOM given 
by Eqs. (j!])-©. All the curves were constructed for the 
RRW06 model and for some values of fi. Notice that 
the dependence of p s and $ with temperature becomes 
smoother when Gy increases. This is consistent with 
our findings that the vector interaction in PNJLv models 
favors the crossover in the quark phase diagram (see Fig. 
[T]). Note also that, the increase of Gy decreases the quark 
density for fixed [i and T values, as expected from the 
relation p = (/i — p)/2Gy. 

The approaches presented here for the construction of 
the hadron-quark phase transitions, can be understood 
from a qualitatively analysis of the structure of the PNJL 
model. It is known that the differences between the NJL 
and the PNJL models, are the modification in the Fermi- 
Dirac distributions and the inclusion of a Polyakov poten- 
tial, W(<I>, $*, T) , in the equations of state. Regarding 
the statistical distributions, the extreme case in which 
$ = $* = leads, as one can verify from Eq. (J3J for 
Gy = 0, to F(k,T,fJ,,^ = $* = 0) = [ e 3(£-M)/T + !]-!_ 



4 



The same result is obtained for the anti-quarks distribu- tions by replacing fi by — fj,. 




Notice that in the case of total confinement, i.e. $ = 
$* = 0, the PNJL model gives rise to the Fermi-Dirac 
statistics of a three-quark cluster, which can be seen 
as a prototype of hadrons. The PNJL model embod- 
ies the formation of hadronic degrees of freedom. From 
this point of view, the interaction has to recognize the 
hadronic formation, and should account for a more realis- 
tic description of the case in which the quarks are totally 
confined. Therefore, the interaction between the clusters 
should also contain the characteristic repulsion present in 
the force between hadrons. This feature is qualitatively 
built in PNJL models, by including the term Gy()/'7 (I ^) 2 
in the Lagrangian density (PNJLv model). 

In the RMF-PNJL approach of the hadron-quark phase 
transition with no vector interaction in the PNJL model, 
all the repulsion is restricted to the hadronic sector. 
What our results reveal is that the effect of such re- 
pulsion in the hadron-quark transition, can be repro- 
duced in PNJLv models by adjusting the vector inter- 
action strength. However, we should point out that at 
high (i and low T, in the region where the baryonic de- 
grees of freedom dominate, models formulated only with 
quark degrees of freedom are not strictly valid (see e.g. 
[37j|). The PNJLv model should also be limited in the 
same sense, while it still incorporates the nuclear repul- 
sion close to the hadron-quark phase transition. As it is 
known (Iol | , by increasing the vector interaction strength 
in the PNJLv model the critical end-point (CEP) disap- 
pears, while in the RMF-PNJL model in all cases it still 
remains. We stress that the vector interaction strength, 
in our model, was fitted in a region of T x fj, far from the 



RMF-PNJL critical end-point (see Fig. [T]). The present 
version of the PNJLv model should be improved to ac- 
count for a possible CEP. For instance, by including den- 
sity and T dependence on the model parameters, in order 
to incorporate more effects from QCD. In particular, the 
recent calculation of the Polyakov quark-meson model 
with the functional renormalization group method [38j . 
already shown that the /^-dependence of To affects the 
QCD phase diagram at finite fx. 

As a last remark, we point out that the construction of 
the hadron-quark phase transition is still an open ques- 
tion, from the experimental point of view. Experiments 
planned to occur in the new facilities, such as the Facil- 
ity for Antiproton and Ion Research (FAIR) [39[ at GSI, 
and the Nuclotron-bascd Ion Collider Facility (NICA) 
[4tj| at the Joint Institute for Nuclear Research (JINR), 
will be performed to reach the highly compressed matter 
covering thus, as far is possible, the strongly interacting 
matter phase diagram. For this reason, it is always im- 
portant to present and discuss the different predictions 
of effective models in the intermediate T and \i values of 
the hadron-quark phase diagram, where the lattice QCD 
calculations still can not reach. 
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